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Driven quantum systems, described by Hamiltonian Ti.{Q, P, x{t)) where x{t) is a time dependent 
parameter, are of interest in mesoscopic physics (quantum dots), as well as in nuclear, atomic 
and molecular physics. Such systems tend to absorb energy. This irreversible effect is known 
as dissipation. More generally, x may be a dynamical variable, where the total Hamiltonian is 
T-Lo{x,p) + Tl{Q,P;x). In such case the interaction of {x,p) with the environmental degrees of 
freedom (Q, P) leads to dephasing as well as to dissipation. It should be emphasized that even 
few {Q, P) degrees of freedom can serve as a miniature heat bath, provided they have chaotic 
dynamics. We shall introduce a general framework for the analysis of dissipation and dephasing, 
and we shall clarify the tight connection to recent studies of quantum irreversibility (also referred 
to as "Loschmidt echo" or as the "fidelity" of quantum computation). Specific model systems that 
will be presented are: particle in a box driven by moving a wall, and particle in a box/ring driven by 
electro-motive- force. These two examples are related to studies of nuclear friction and mesoscopic 
conductance. Specific issues to be discussed are the limitations of kinetic theory, the capabilities of 
linear response theory, and the manifestation of non-perturbative quantum-mechanical effects. In 
particular we shall explain that random matrix theory and the semiclassical theory lead to different 
non-perturbative limits. 
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I. INTRODUCTION 

In the following lectures we are interested in systems 
that are described by a Hamiltonian 7i((5, P', x{t)), where 
{Q,P) is a set of (few) canonical coordinates, and x{t) 
is a time dependent parameter. We further assume that 
in the time independent case {x{t) = const) the classical 
motion is chaotic. 

The Quantum mechanical (QM) study of classically 
chaotic systems is known in the literature as "quantum 
chaos". The theory of time independent Hamiltonians, 
in particular their spectral properties, is well documented 

^1 «2J ■ -But the QM theory of driven chaotic systems is 
still a missing chapter. 

To avoid misunderstanding we refer here to systems 
that are chaotic also in the absence of driving. By "driv- 
ing" we mean that we assume Hamiltonian ^{Q, P; x{t)) 
where x{t) is time dependent, rather than constant. Un- 
like the case of time independent systems, in case of 
driven systems the energy distribution evolves with time. 
What one needs is a theory for this evolution. We shall 
see that various notions, such as "dissipation", "irre- 
versibility" and "fidelity" just emphasize particular as- 
pects of this evolution. 

The available theory for driven systems in the quan- 
tum mechanical literature is based mainly upon first or- 
der perturbation theory, supplemented by higher order 
(sometimes non-perturbative) corrections. Depending on 
the "rate" of the driving, one commonly distinguishes 
between the "adiabatic" and the "Fermi golden rule" 
regimes 0, |^. The latter is known also as the "linear 
response theory" (LRT) or as the "Kubo-Greenwood" 
regime. 

There are circumstances where first order perturbation 
theory cannot serve as a starting point for the analysis of 



a driven system. However the well known examples refer 
to systems that are not chaotic in the absence of driv- 
ing. This includes in particular one-degree-of-freedom 
systems such as the quantum kicked rotator j^] . Our in- 
terest is in generic chaotic systems, therefore we have to 
consider systems that have at least two degrees of free- 
dom. 

The possibility to present a general QM theory for 
driven chaotic systems 0, 0, 0| follows from the 
simple fact that "chaos" leads to universality. This uni- 
versality is captured, to some extend, by random matrix 
theory (RMT). On the other hand we have semiclassi- 
cal methods. We shall see that there is a clash between 
RMT and semiclassics, and that they lead to different 
non-perturbative limits jl3lll4|. 



II. MODEL SYSTEMS 

The classic example for a driven system is the piston 
model (Fig.l), where a gas in confined inside a cylinder, 
and X is the position of the piston. Our interest is in the 
case where we have "one particle gas". [Note however 
that if we know to solve the problem for one particle, 
then automatically we can get the solution for many non- 
interacting particles] . 

The ID-box version of this model (Fig. 2a) is known in 
the literature as the " infinite- well" problem with moving 
wall [TH . Some limited aspects of this problem have 
been discussed in the literature in connection with the 
Fermi acceleration problem ^3 ■ 

A 2D-box variation of the "piston model" is presented 
in Fig. 2b. Here we have stadium shaped billiard, and the 
the parameter x controls the deformation of the bound- 
ary. Two other variations of the same model are pre- 
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sented in Fig. 3, where the box has the shape of a gener- 
ahzed Sinai biUiard. 

In the examples so far the parameter x controls the 
shape of the " box" , and V = x has the interpretation of 
wall velocity. The interest in such systems has emerged 
long time a go in studies of nuclear friction (one-body 
dissipation) |l8l [T9| . A renewed interest is anticipated in 
mesoscopic physics where the shape of a quantum dots 
can be controlled by gate voltages. [Note that in the 
nuclear physics context the shape is close to spherical, 
while a quantum dot is typically strongly chaotic] . 

We can create driving by changing any parameter 
(or field). In Fig. 4 the driving is achieved by chang- 
ing the perpendicular magnetic field. Fig. 4a assumes 
"quantum dot geometry" with homogeneous magnetic 
field, while Fig. 4b assumes Aharonov-Bohm ring geom- 
etry with magnetic flux that goes via the hole. Let us 
define x as the total magnetic flux. In such case V — x 
is the electro-motive force (measured in Volts) which is 
induced in the ring according to Faraday law. 

If the variations of the parameter x are classically 
small, then we can linearize the Hamiltonian as follows 

■H{Q,P;x)^nQ{Q,P)+xW{Q,P) (I) 

where without loss of generality we have assumed that 
X = is the typical value of x. For generic systems 
(which means having smooth Hamiltonian that generates 
a classically chaotic motion), the representation of W, 
in the ordered Tio determined basis, is known to be a 
banded matrix (for details see the next section) . A simple 
example can be found in |20l | where 

H(Q, P; x) = i(p2+p| + Q\+Ql) + (1 + x) • QlQl (2) 

This Hamiltonian describes a particle moving inside a two 
dimensional anharmonic well (2DW). The shape of the 
2DW in controlled by the parameter x. The perturbation 
is yV((5, P) = QiQii 8.nd its matrix representation B„m 
is visualized in the inset of Fig. 7. 

The above discussion of generic Hamiltonian models, 
such as the 2DW model, motivates the definition of a 
simple artificial model Hamiltonian, that has the same 
characteristics: This is Wigner model pH |2^ . In the 
following definition of Wigner model we follow closely the 
notations of p^. In the standard representation Tig = 
E is a diagonal matrix whose elements are the ordered 
energies with mean level spacing A, and W = B 

is a random banded matrix with non- vanishing couplings 
within the band < |n— m| < b. These coupling elements 
are zero on the average, and they are characterized by the 
variance a = {{\Rnm\'^)Y'^ ■ Hence the Hamiltonian is 

7^ = E + a;B (3) 

This artificial model can serve as a reference case for 
testing various theoretical ideas. Moreover, it has been 
conjectured that such model captures some generic fea- 
tures of chaotic systems. [Note that most of the RMT 
literature deals with simplified versions of Wigner model, 
where the bandwidth equals to the matrix size]. 



III. QUANTUM CHAOS 

The notion of chaos in classical mechanics implies that 
few degree of freedom system, such as the Sinai billiard 
system, exhibit stochastic-like behavior. This is in con- 
trast to the out-of-date idea that stochasticity and irre- 
versibility are the outcomes of having (infinitely) many 
degrees of freedom. Chaos means that the motion (eg 
Fig. 5) has exponential sensitivity to any perturbation or 
change in initial conditions. Another way to character- 
ize a chaotic motion is by its continuous power spectrum 
(see Fig. 6). This should be contrasted with integrable 
motion which is characterized by a discrete (rather than 
continuous) set of frequencies. 

For sake of later analysis it is useful to define the 
"power spectrum" of the motion specifically as follows. 
Let {Q{t), P{t)) be an ergodic trajectory that is gener- 
ated by the time independent Hamiltonian Ti.(Q, P; x). 
We can define a fluctuating quantity JF(i) — —dH/dx. 
In case that x is the displacements of a wall element 
(eg Fig. 3b), the fluctuating T{t) has the meaning of 
" Newtonian force" . In case that x is the magnetic flux 
(eg Fig. 4b), the fluctuating T{t) has the meaning of 
"electric current". In case of the 2DW model we get 
jc-(t) = -W{Q{t),P{t)) = -Qi{tfQ2{tf. The correla- 
tion function of the fluctuating T{t) will be denoted by 
C (r) and the power spectrum of the fluctuations will be 
denoted by C{ijj). The latter is the Fourier transform of 
the former. The variance of the fluctuation is C(0), the 
intensity of the fluctuations is defined as C(0), and the 
correlation time is denoted by Td- 

It is clear that upon quantization we no longer have 
chaos. Still, the question arise what are the fingerprints 
of the classical chaos on both the spectral properties of 
the system, and also on the structure of the eigenstates. 
This problem was the focus of intensive studies during the 
last decade Q, , and it has important applications in 
mesoscopic physics |3 12^ . 

An important observation of " quantum chaos" studies 
is that Quantum Mechanics introduce two additional en- 
ergy scales into the problem (rather than only one). We 
can take the 2DW model as a generic example. After 
rescaling of the classical parameters of the model, we are 
left with one dimensionless parameter (the dimensionless 
energy). This parameter controls the nature of the classi- 
cal dynamics. Upon quantization we have two additional 
(dimensionless) parameters. One energy scale is obvi- 
ously the mean level spacing A, which is proportional to 
fi^ . The other energy scale is A;, = Ti/tcIj where t^i is the 
classical correlation time that characterizes the (chaotic) 
dynamics. If h is small then the two energy scales are 
very different (A ^ A;,). 

The significance of the energy scale Ah is a central is- 
sue in " quantum chaos" . It turns out that the statistical 
properties of the energy spectrum are universal on the 
sub-?i scale, and obey the predictions of RMT. On the 
other hand, on large energy scale (compared with A^), 
non-universal (system specific) features manifest them- 
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selves |2g. These features are the fingerprints of the 
underlying classical dynamics. In the context of ballistic 
quantum dots, which are in fact billiard systems, Af, is 
also know as the " Thouless energy" . 

There is another way in which the energy scale Ab 
manifests itself. Let W{Q, P) be some observable, and 
consider its matrix representation B„„i in the basis which 
is determined by the (chaotic) Hamiltonian. An exam- 
ple is presented in Fig. 7. It can be argued |23| that 
B„„i is a banded matrix, and that the bandwidth is 
b = Af,/A. This is based on a remarkably robust semi- 
classical expression that relates the bandprofile to the 
classical power spectrum: 

We can apply this semiclassical relation to the case where 
yV is the "perturbation" as defined in Eq.Q. This leads 
to the interpretation of A^ as the largest "distance" in 
energy space that can be realized in a first-order tran- 
sition. We can also use the semiclassical relation in 
reverse, in order to find/define the classical correlation 
function that corresponds to a quantum-mechanical ma- 
trix Hamiltonian. In case of the standard Wigner model 
we get C(t) — b(T'^sm{T / Tci) / {t / Tci) , with the correlation 
time Tci = h/ (6A). 



IV. PARAMETRIC EVOLUTION 

A more recent development was to consider a para- 
metric set of Hamiltonians, namely T-l{Q, P;x) where x 
is a parameter as in the examples of Section 2. For each 
value of X we can diagonalize the Hamiltonian, leading to 
set of (ordered) eigen-energies E„ (x) , as in the schematic 
illustration of Fig. 8. The corresponding eigenstates will 
be denoted by \n{x)). Their parametric evolution can be 
characterized by the parametric kernel 

P{n\m) — \ {7i{x)\m{xQ))\^ (5) 

We shall use the notation P(r) — P{n ~ m) — P(n|m), 
with implicit average over the reference state m. We shall 
refer to P{r) as the "average spreading profile". This is 
in fact, up to scaling, the LDOS (local density of states, 
also known as strength function). 

Let us characterize the perturbation by the quantity 
5x = x — xq. The interesting question is how P{r) evolves 
as we increase the perturbation 5x. For the Wigner 
model the answer is known long ago [^|2^|2^|23. P(r) 
has a standard perturbative structure for very small 5x. 
For larger 5x it becomes a chopped Lorentzian, and for 
even larger 5x it becomes a semicircle. We shall denote 
the border between the standard perturbative regime 
and the Wigner regime by 5xc, and the border between 
the Wigner regime and the non perturbative (semicircle) 
regime will be denoted by (5xp,t . The explicit expressions 



are: 

5x, ^ % oc (6) 

where d is the number of freedoms {d — 2 for billiards) . 
In order to determine the fi dependence we have used 
the semiclassical relation Eq.Q, and the proportionality 
A cx ?i . Note that the latter relation, known as Wyle 
law, is significant for the determination of Sx^- In con- 
trast to that feprt cx ?i is in fact independent of A. 

The generalization of Wigner scenario has been the 
subject of our recent research [23, 113, l^- In the general 
case the standard perturbative structure evolves into a 
" core-tail structure" , while for large 5x it becomes purely 
non-perturbative. In the standard perturbative regime 
{5x <C 5xc) most of the probability P{r) is concentrated 
in one level (n = to). In the extended perturbative 
regime most of the probability is concentrated within a 
"core" whose width is typically V{5x) ^ {5x/5xc)^. The 
"core" is the non-perturbative component which arise 
due to non-perturbative mixing of nearby levels. The 
"tail" is the outer perturbative component which is cre- 
ated by first order transitions. 

The extended perturbative regime is defined by the re- 
quirement of having separation of energy scales T{5x) <C 
Aft. This condition is trivially satisfied in the "stan- 
dard perturbative regime" where F ^ A. The condition 
r(i5x) <C Ah is violated in the non-perturbative regime 
{5x 3> (^Xprt), which in fact leads to the determination 
of (5xprt as in Eq.©. The theory for P{r) in the non- 
perturbative regime is not complete yet. However, it can 
be argued js^l that if 5x is large enough, then P(r) be- 
comes of semiclassical nature [sj • The case of billiards 
with shape deformation re quir es special considerations 
and is of particular interest [S^, ISJ . 

It is important to realize that the border of the stan- 
dard perturbative regime {5xc) is related to the energy 
scale A, while the border of the extended perturbative 
regime (Jxp^t), which leads to the identification of the 
non-perturbative regime, is related to the bandwidth A;,. 

V. TEMPORAL EVOLUTION 

After considering the parametric evolution, the next 
logical stage is to consider the actual (temporal) evolu- 
tion which is generated by the time dependent Hamilto- 
nian Ti.{Q , P, x{t)) . Then, in complete analogy, we can 
ask how the energy scales A and A;, are reflected in the 
actual evolution. We postpone the discussion of the lat- 
ter issue to Section 11. 

The purpose of the present and next sections is to de- 
fine what does it mean " driving" , and how do we quantify 
the temporal evolution. Without loss of generality we as- 
sume x{0) = 0. We would like to consider the following 
driving schemes: 
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• Linear driving 

• One pulse driving cycle 

• Periodic driving 

• Driving reversal scenario 

• Time reversal scenario 

In the next section we define the various schemes, some of 
which are also illustrated in Fig. 9. The evolution is char- 
acterized by the obvious generalization of Eq.jSJ, namely 

Ptin\m) = \{n{x{t))\U\mix{0)))\^ (8) 

Here U is the evolution operator, with implicit depen- 
dence on the time t. The parametric kernel Eq.|SJl can 
be regarded as corresponding to the " sudden" limit where 
[/ « 1. As in the parametric case we can define an aver- 
age spreading profile -Pt(r), where r = n ~ m. 

There are various practical possibilities available for 
the characterization of the distribution Pt{r). It turns 
out that the major features of this distribution are cap- 
tured by the following three measures: 

V{t) = Pt(r = 0) (9) 
r{t) = 50% probabihty width (10) 



SE{t) = 




The first measure is the survival probability Vit). The 
second measure T(t) is the energy width of the central r 
region that contains 50% of the probability. [For simplic- 
ity of presentation we use here a loose definition of F as 
far as prefactors of order unity are concerned.] Finally 
the energy spreading 5E{t) is defined as the square-root 
of the variance. 

It is important to realize that the above three measures 
give different type of information about the nature of the 
energy spreading profile. In particular the indication for 
having a core-tail structure is: 

T{t) < 5E{t) < Ab (12) 

The core-tail structure (eg chopped Lorentzian) is charac- 
terized by a " tail" component that contains a vanishingly 
small probability but still dominates the variance. [Note 
that Af, = oo, as in the case of un-chopped Lorentzian, 
would imply SE{t) — oo irrespective of F.] In contrast to 
that a typical semiclassical spreading profile (as well as 
Wigner's Semicircle) is characterized by 

6E{t) - r{t) > Af, (13) 

In the latter case, in order to avoid confusion, it is bet- 
ter not to use the notation T{t). [The notation F has 
been adopted in the common diagrammatic formulation 
of perturbation theory. This formulation is useful in the 
extended perturbative regime in order to derive Wigner's 
chopped Lorentzian. In the non-perturbative regime this 
formulation becomes useless, and therefore the signifi- 
cance of F is lost.] 



VI. DRIVING SCHEMES 

Linear driving means x(t) — Vt, where is a con- 
stant. In such a scenario obviously x{t) ^ x{0). Still it 
is convenient to assume that the chaotic nature of the 
dynamics is independent of x, and that changes in x are 
not associated with changes in phase space volume (no 
conservative work is being done). This is manifestly the 
case for the systems which are illustrated in Fig. 3b and 
Fig. 4b. [Note that the standard Wigner model does not 
have x invariance property, and therefore the analysis of 
linear driving for the Wigner model is an ill defined prob- 
lem. Attempts to overcome this difficulty lead to certain 
subtleties (llij]. 

For all the other driving schemes we assume, without 
loss of generality, that a;(0) = x{T) = 0, where T is the 
period of the driving cycle. The simplest driving scheme 
is a rectangular pulse x{0 < t < T) = A, which is charac- 
terized by its amplitude A. Does the study of rectangular 
pulses constitute a good bridge for developing a general 
theory for driven systems? The answer is definitely not. 
An essential ingredient in the theory of driven system is 
the rate V in which the parameter x is being changed in 
time. Therefore, it is important to consider, for exam- 
ple, a triangular pulse (Fig. 9b). Such pulse is character- 
ized by both amplitude A and driving rate V — 2A/T. 
More generally one may consider (Fig. 9c) a train of such 
pulses (= periodic driving). In particular one may con- 
sider the usual sinusoidal driving x{t) = Asm{nt) where 
il = 2tt/T. In the later case we can define the root- mean- 
square driving rate as V = flA/^/2. 

In all these cases we ask, in complete analogy with 
the parametric case, what is the evolution of the energy 
distribution Pt{f). Now the evolution is with respect 
to time, rather than with respect to 5x. The different 
scenarios are distinguished by the choice of U . We shall 
use the notation f7[a;yi] in order to denote the evolution 
operator that corresponds to driving scheme x — XA{t). 

The case of rectangular pulse is known in the liter- 
ature as "wavepacket dynamics" [s^]- The particle is 
prepared in an eigenstate of the T-Lq Hamiltonian, while 
the evolution is generated by the perturbed Hamiltonian 
H. = H{Q,P;x = A). We may consider more compli- 
cated schemes of pulses. For example rectangular pulse 
+A followed by another rectangular pulse —A. The ques- 
tion here is whether the second pulse can compensate 
the effect of the first pulse. We call such scheme "driv- 
ing reversal" . The evolution operator can be written as 
U — L/[a;(rev)]J7[a;] where x = A represents the rectangu- 
lar pulse, while a::(rev) = —A is the reversed pulse. The 
case of triangular pulse can be regarded as another par- 
ticular variation of driving reversal. In the latter case x 
represents linear driving with velocity Y , and x(rev) is 
the reversed driving process with velocity —V. 

The case of " driving reversal" should be distinguished 
from "time reversal" scheme. The latter notion is ex- 
plained in the next section. 
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VII. TWO NOTIONS OF IRREVERSIBILITY 



rectangular pulse. The survival amplitude is defined as 



There are two distinct notions of irreversibility in sta- 
tistical and quantum mechanics. One is based on the 
"piston model" paradigm (PMP), while the other [s^l is 
based on the "ice cube in cup of hot water" paradigm 
(ICP). The recent interest [sl IM El S IM 113 in 
"quantum irreversibility" is motivated by its relevance 
to quantum computing. 

In the PMP case we say that a process is reversible if 
it is possible to "undo" it by driving reversal. Consider 
a gas inside a cylinder with a piston (Fig.l). Let us 
shift the piston inside. Due to the compression the gas is 
heated up. Can we undo the "heating" simply by shifting 
the piston outside, back to its original position? If the 
answer is yes, as in the case of strictly adiabatic process, 
then we say that the process is reversible. 

In the ICP case we consider the melting process of an 
ice cube. Let us assume that after some time we reverse 
the velocities of all the molecules. If the external condi- 
tions are kept strictly the same, we expect the ice-cube to 
re-emerge out of the water. In practice the external con- 
ditions (fields) are not exactly the same, and as a result 
we have what looks like irreversibility. 

The mathematical object that should be considered in 
order to study PMP is just Pt{r) for a driving scheme 
that involves " driving reversal" . Namely, as discussed in 
the previous section, the evolution operator is 



F{t) = (*o|C/[A]|^'o) 



(16) 



U 



U[x{rev)] U[x] 



(14) 



The mathematical object that should be considered in or- 
der to study ICP is again Pt (r) , but with driving scheme 
that involves " time reversal" . Namely, the evolution op- 
erator is defined as 



U = U[xb]-^ U[xa] 



(15) 



In the latter case, if the external conditions are in full 
control {xb — xa), then we have complete reversibility 

It is also important to define precisely what is the mea- 
sure for quantum irreversibility. This is related to the 
various possibilities which are available for the charac- 
terization of the distribution Ft (r) . The prevailing possi- 
bility is to take the survival probability V{t) as a measure 
jsJl- Another possibility is to take the energy spreading 
6E{t) as a measure. The latter definition goes well with 
the PMP, and it has a well defined classical limit. Irre- 
versibility in this latter sense implies diffusion in energy 
space, which is the reason for having energy absorption 
(dissipation) in driven mesoscopic systems (see section 9). 



VIII. WAVEPACKET DYNAMICS, SURVIVAL 
PROBABILITY AND FIDELITY 



exp 



The survival probability is 'P{t) — \F{t)\'^, in consistency 
with the definition of Eq.©. F{t) is manifestly a Fourier 
transform of the LDOS, and therefore we can immedi- 
ately draw a conclusion js^l that the nature of the dy- 
namics is different depending on the parametric regime to 
which the amplitude A belongs. If P{r) have a core-tail 
Lorentzian structure, then we get an exponential decay 
Vit) — exp(— Pi). On the other hand if P{r) has a semi- 
classical structure, then the decay ofV{t) is non- universal 
(system specific). 

A similar picture arise in recent studies of the survival 
probability for "time reversal" driving scheme. Here one 
defines the fidelity amplitude as 



F{t) - {^o\U[A]'^U[0]\^o) 



(17) 



The fidelity, also known as Loschmidt echo, is defined as 
Vit) — \F{t)\'^. The situation here is more complicated 
compared with Ea. (|16|l because we have two LDOS func- 
tions [s^: one is the Tio weighted LDOS, and the other 
is the ^'o weighted LDOS. The two LDOS functions co- 
incide only if ^0 is an eigenstate of Hq. In the latter 
case the F{t) of Ea. (|17|) reduces (up to phase factor) to 
Eq. (|16|) . It turns out that in case of Eq. (|17l) there is no 
simple Fourier Transform relation between F{t) and the 
LDOS functions. However, the picture "in large" is the 
same as in the case of Ea. (|16|) js^- Namely, one has 
to distinguish between three regimes of behavior: In the 
standard perturbative regime (A < Sxc) one typically 
encounters a Gaussian decay (3^]; In the Wigner regime 
(also called FGR regime) one typically finds Exponen- 
tial decay [s^ : And in the non-perturbative regime one 
observes a semiclassical perturbation-independent "Lya- 
punov decay" [ssf. 

The study of the survival probability, as described 
above, is only one limited aspect of the temporal evolu- 
tion. The more general object that should be considered 
is Pt{r) as defined in section 5. The major features of this 
time evolution are captured by the three measures that 
we have defined in Ea.l|!?|l- lll|) . In Fig. 10 we display nu- 
merical simulations of wavepackct dynamics for the 2DW 
model. The energy spreading SE{t) is plotted as a func- 
tion of time. The first panel is the classical simulation, 
which in fact coincides with the "linear response" calcu- 
lation. The input for the LRT calculation is C(t), and 
the result is proportional to the amplitude A. Namely, 



SE{t) = A X ^/2{C{0)~C{t)) 



(18) 



Driving schemes with rectangular pulses are the sim- 
plest for both analytical and numerical studies. It is eas- 
iest to consider the survival probability in case of a single 



In the second panel we display the results of the quan- 
tum mechanical simulations. For smaller h values the 
agreement with the classical LRT calculation is better. 
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Finally, in the third panel we repeat the quantum me- 
chanical simulations with a sign randomized Hamilto- 
nian. This means that we take Eq. Q , and we randomize 
the sign of the ofF-diagonal terms. The bandprofile, and 
hence C(w) are not affected by this procedure, which 
implies that the LRT calculation gives exactly the same 
result. But now we see that as h becomes smaller the 
correspondence with the classical result becomes worse. 
Specifically: In (a) and (b) we see a crossover from bal- 
listic spreading {5E oc t) to saturation {SE ~ const) as 
implied by Ea. H18|) . Only one time scale (r^i ~ 1) is in- 
volved. In (c), in contrast to that, we see that as ?i ^ 
an intermediate stage of diffusion {SE cx Vt) develops. 

How can we explain the above results. Obviously we 
see that for small h we cannot trust LRT. What in fact 
happens is that we have a crossover from the perturbative 
regime {A < 6xp^t{h)) to the non-perturbative regime 
{A > Sxp^t{h)). In the latter case we get either semiclas- 
sical behavior, or RMT behavior. In other words, random 
matrix theory and the semiclassical theory lead to differ- 
ent non-perturbative limits. In the semiclassical case the 
crossover from LRT behavior to non-perturbative behav- 
ior cannot be detected by looking on SE{t). Still the 
crossover can be detected by looking on T{t). See 'T4t| for 
details. 



IX. DIFFUSION IN ENERGY SPACE AND 
DISSIPATION 

In the following sections we discuss the case of either 
linear or periodic driving. In such case the long time be- 
havior of the system is characterized by diffusion in en- 
ergy space. Associated with this diffusion is a systematic 
increase of the average energy. This irreversible process 
of energy absorption is known as " dissipation" . 

There is a satisfactory classical theory for dissipation 
[4l|. Some of the mathematical details are subtle, but 
the overall physical picture is quite simple. Without loss 
of generality the main idea can be explained by referring 
to the bilhard example of Fig. la. The particle executes 
chaotic motion, and we may say that each collision has 
roughly equal probability to be with either the inward- 
going or with the outward-going wall. As a result the 
particle either gain or loose kinetic energy. Thus, the 
dynamics in energy space is like random walk, and it can 
be described by a diffusion equation. Thus we see that 
due to the chaos we have stochastic-like energy spreading. 

This classical diffusion process is irreversible in the 
PMP sense. Let us assume that we start with a mi- 
crocanonical distribution that has definite energy. If, af- 
ter some time, we reverse the velocity of the walls, then 
obviously we do not get back the initial microcanonical 
distribution. 

The effect of dissipation is related to the irreversible 
stochastic-like diffusion in energy space. If the diffusion 
rate were the same irrespective of the energy, then ob- 
viously the average energy would be constant. But this 



is not the case. The diffusion is stronger as we go up 
in energy, and as a result the diffusion process is biased. 
Thus the average energy systematically grows with time, 
and one can derive a general diffusion dissipation relation 
Q: 

where g{E) is the density of states, and p{E) is the proba- 
bility distribution (eg microcanonical, canonical or Fermi 
occupation). The diffusion picture is generally valid in 
the classical case, and it is typically valid also in the 
quantum mechanical case. [The issue of dynamical local- 
ization due to strictly periodic driving is important 
for driven ID system, but not so important in the case 
of driven chaotic systems |Tll|.] 

There is a simple linear response (Kubo) expression, 
that relates the diffusion coefficient to the power spec- 
trum C{u>) of the fluctuations: 

= ^Cin) X (20) 

The diffusion law for short times is 5E{t) — \/2lJ^. This 
expression is completely analogous to Ea. H18|) . In both 
cases the spreading is proportional to the amplitude A. 
[Recall that for periodic driving we define V = QA/V^. 
In the special case of linear driving the spreading is pro- 
portional to v.] Moreover, as in the case of wavepacket 
dynamics, the LRT result is the same classically and 
quantum-mechanically. But again, as in the case of 
the wavepacket dynamics, the validity regime of LRT in 
the quantum mechanical case is much smaller (see sec- 
tion. 11). 

If we combine the above Kubo expression with the 
diffusion-dissipation relation we get 

l{n)^^,in)xv' (21) 

where /i is related to the power spectrum of the fiuctu- 
ations. Thus we get a fluctuations-dissipation relation 
|43 |. The standard "thermal" fluctuation-dissipation re- 
lation ^(0) = C{0)/{2kBT) is obtained from Eq.^ in 
case of canonical p{E). 

Standard textbook formulations 01 takes linear re- 
sponse theory together with thermal statistical assump- 
tions as a package deal. Our presentation provides a 
more powerful picture. On the one hand we can discuss 
non-equilibrium situation using LRT combined with an 
appropriate version of the diffusion-dissipation relation. 
On the other hand, we may consider situation where LRT 
does not apply. In such case we may get some (non- 
perturbative) result for the diffusion, and later use the 
diffusion-dissipation relation in order to calculate the dis- 
sipation rate. 
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X. BEYOND KINETIC THEORY 

The coefficient fi in Ea. (|21|l is called the "dissipation 
coefficient". In the case where x is the displacement of 
a wall element, it is also known as "friction coefficient", 
and in the case where a; is a magnetic flux it can be called 
" conductance" . 

Having dissipation rate proportional to is known as 
"ohmic" behavior. In case of "friction" it is just equiva- 
lent to saying that there is a friction force proportional 
to the velocity V, against which the wall is doing me- 
chanical work. This mechanical work is " dissipated" and 
the gas is " heated up" in a rate proportional to V'^ . 

In case of "conductance" we may say that there is a 
drift current proportional to the voltage V. This is in fact 
"Ohm law". The dissipated energy can either be accu- 
mulated by the electrons (as kinetic energy), or it may be 
eventually transfered to the lattice vibrations (phonons). 
In the latter case we say that the ring is " heated up" . The 
rate of the heating goes like which is in fact "Joule 
law". 

The traditional approach to calculate /i is to use a ki- 
netic picture (Boltzmann) which is based on statistical 
assumptions. This leads in case of friction to the "wall 
formula" 

N 

Hin) = —mv^A„^n, (22) 

where N is the number of gas particles (let us say iV = 1), 
and Ve = ^/2E/m. We also use the notations V^ox for 
the volume of the box, and A„^ii, for the effective area of 
the mo ving walls. In the latter we absorb some geometric 
factors [ij . Application of the traditional kinetic (Boltz- 
mann) approach in case of conductance leads to "Drude 
formula" : 

where A^^t is the area of the " quantum dot" , and t^^i is 
the average time between collisions with the walls. 

Using Linear response theory (Kubo formula), as de- 
scribed in the previous section, we can go beyond the 
over-simplified picture of kinetic theory. That means to 
go beyond Boltzmann picture. Below we explain under 
what assumptions we get the "traditional" kinetic ex- 
pressions, and what in fact can go wrong with these as- 
sumptions. 

The interest in friction calculation has started in stud- 
ies of "one body dissipation" in nuclear physics 0,^3- 
The "wall formula" assumes that the collisions are totally 
uncorrelated. In such case the power spectrum C{uj) of 
!F{t) is like that of white noise (namely "flat"). By in- 
spection of Fig. 6 we can see that this assumption is ap- 
parently reasonable in the limit of very strong chaos. But 
it is definitely a bad approximation in case of weak chaos. 
The dynamics of chaotic system is typically characterized 
by some dominant frequencies. Therefore we have rel- 
atively strong response whenever the driving frequency 



matches a "natural" frequency of the system. This can 
be regarded as a classical (broad) resonance. By inspec- 
tion of Fig. 6 we see that a particular feature is having 
such resonance around uj — 0. This type of resonance, 
due to bouncing behavior, can be regarded as a "classical 
diabatic effect" 

Even if the chaos is very strong, the "white noise" as- 
sumption is not necessarily correct: In .45,, 47] we explain 
that for special class of deformations (including transla- 
tions, rotations and dilations) the low frequency response 
is suppressed, irrespective of the chaoticity. This is illus- 
trated numerically in Fig. 11. 

In case of Drude formula the fluctuating ^{1) has the 
meaning of " electric current" , and therefore the power 
spectrum C(u;) is the Fourier transform of the current- 
current (or one may say velocity-velocity) correlation 
function. Assuming that the velocity-velocity correla- 
tion function decays exponentially in time, one obtains 
the Drude result. A careful analysis of this assumption, 
and its relation to the "white noise approximation" of 
the " wall formula" , can be found in Section 6 of . 
Fig. 12 displays a numerical example. We clearly see non- 
universal deviations from the Drude expression, which 
reflect the specific geometry of the " quantum dot" . 



XI. NON-PERTURBATIVE RESPONSE 

In the classical case, assuming idealized system, the 
crossover to stochastic energy spreading involves only one 
time scale, which is t^,. Gaussian-like spreading profile is 
obtained only for time t much larger than r^,,. For short 
times we can use linear analysis in order to calculate the 
spreading profile. However, this analysis has a breaktime 
that we call tf„(T^), where V is the rate in which x 
is being changed. For long times {t ^ r^i) we can use 
stochastic picture. Classical LRT calculation of the dif- 
fusion is valid only if the crossover to stochastic behavior 
is captured by the short time analysis. This leads to the 
classical slowness condition r^i ^ t!,^{V) which we as- 
sume from now on. See specific examples in Sections 13 
and 14. 

In the quantum mechanical case we follow a similar 
reasoning. The linear analysis is carried out using pertur- 
bation theory. We have presented a careful analysis 
to determine the breaktime tprt(T^;?i) of this analysis. It 
turns out that this breaktime is not related to the mean 
level spacing A, but rather to the much larger energy 
scale Af). 

In complete analogy with the classical analysis, it turns 
out that the validity of LRT in the quantum domain is 
restricted by the condition <C tp^t{V;h). If this in- 
equality is not satisfied, then we say that we are in a 
non-perturbative regime. It is important to realize that 
the h —^ limit is a non-perturbative limit. This means 
that the semiclassical regime is contained within the non- 
perturbative regime. 

In the simple examples that are discussed in Sections 
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13 and 14, the non-perturbative regime is in fact a semi- 
classical re gim e. This coincidence does not hold in gen- 
eral p. IT(iTl3|. In case of RMT models, obviosly we do 
not have a semiclassical limit. In such models the non- 
perturbative response deviates significantly from Kubo 
formula (Fig. 13). The interest in such models can be 
physically motivated by considering transport in quan- 
tized disordered systems. Whether similar deviations 
from Kubo formula can be found in case of quantized 
chaotic systems is still an open question [TH . In any 
case, it is important to remember that the rate of dissi- 
pation is just one aspect of the energy spreading process. 
Even if Kubo formula does not fail (thanks to quantum- 
classical correspondence), still there are other features of 
the dynamics that are affected by the crossover from the 
perturbative to the non-perturbative regime. For exam- 
ple: in Sec. 19 we are going to show that different results 
are obtained for the dephasing time, depending whether 
the process is perturbative or non-perturbative. 

We can express the condition for being in the non- 
perturbative regime as 0| 

y » ^ (24) 

Tcl 

The expression in the right hand side scales like h, which 
reflects that this condition is related to and not to A. 
In the next section we discuss the definition of the adi- 
abatic regime (very small V) whose existence is related 
to having finite A. A schematic illustration of the three 
regimes (adiabatic, LRT, non-perturbative) is presented 
in Fig. 14. Some further reasoning 11] allows to define 
the non perturbative regime in case of periodic driving. 
Its location in {A,n) space is also illustrated in Fig. 14. 
Note that for periodic driving we define V = nA/V2. 
The two V — const curves in the (fi, A) diagram repre- 
sent the same conditions as in the case of linear driving. 
Other details of this diagram are discussed in the next 
section, and in [TH . 

XII. ADIABATIC RESPONSE AND QM 
RESONANCES 

Let us assume that we are in the perturbative regime 
(which means that the non-perturbative regime of the 
previous section is excluded) . We ask the following ques- 
tion: can we use the classical Kubo result as an approx- 
imation for the quantum mechanical result? The answer 
is "yes" with the following restrictions: (i) The amplitude 
of the driving should be large enough; (ii) The frequency 
of the driving should be large enough. The two conditions 
arc further discussed below. If they are satisfied we can 
trust the classical result. This follows from the remark- 
able quantal-classical correspondence which is expressed 
by Eq.Q). We have an illustration of this remarkable 
correspondence in Figures 7 and 11. 

Large enough amplitude means A ^ Sxc- One may say 
that large-amplitude driving leads to effective "broaden- 



ing" of the discrete levels, and hence one can treat them 
as if they form a continuum. This is essential in order 
to justify the use of Fermi golden rule (FGR) for a small 
isolated system TTj. Kubo formula can be regarded as 
a consequence of FGR. If the driving amplitude is not 
large enough to "mix" levels, we cannot use FGR, but we 
can still use first order perturbation theory as a starting 
point. Then we find out, as in atomic physics applica- 
tions, that the response of the system is vanishingly small 
unless the driving frequency u; matches energy level spac- 
ing. This is called " QM resonance" . The strips of QM 
resonances are illustrated in the schematic diagram of 
Fig. 14. It is important to realize that higher order of per- 
turbation theory, and possibly non-perturbative correc- 
tions, are essential in order to calculated the non-linear 
response in this regime 0. Still, first order perturbation 
theory is a valid starting point, and therefore we do not 
regard this (non-linear) regime as " non-perturbative" . 

Large enough frequency means uj 3> A/h. The re- 
markable quantal-classical correspondence which is ex- 
pressed by Eq.|01) is valid only on energy scales that are 
large compared with A. If this condition is not satisfied, 
we have to take into account the level spacing statistics 
[isL |4^ . This means that we can have significant dif- 
ference between the quantal LRT calculation, and the 
classical LRT calculation. 

However, this is not the whole story. If V is small 
enough, first order perturbation theory implies " QM adi- 
abaticity". The condition for QM adiabaticity is F ^ 
Sxc/tfi where = 27r?i/A is the Heisenberg time. A 
useful way of writing the QM adiabaticity conditions is: 

" « ^ (^) 

In the adiabatic regime, first order perturbation theory 
implies zero probability to make a transition to other 
levels. Therefore, to the extend that we can trust the 
adiabatic approximation, all the probability remains con- 
centrated in the initial level. Thus, in this regime, as in 
the case of small amplitudes, it is essential to use higher 
orders of perturbation theory, and non-perturbative cor- 
rections (Landau-Zener 0|). Still we emphasize that first 
order perturbation theory is in fact a valid starting point, 
and therefore we do not regard this (non-linear) regime 
as " non-perturbative" . 

XIII. DRIVING BY ELECTRO-MOTIVE FORCE 

Consider a charged particle moving inside a chaotic 
ring. Let x represent a magnetic flux via the ring. If 
we change x in time, then by Faraday law ^ = a; is the 
electro-motive force (measured in Volts). The fluctuating 
quantity J^(t) has the meaning of electric current. The 
variance of the fluctuations is C(0) — (cve/L)'^, where 
We — y^2E/m, and L is the length of the ring. The 
correlation time of these fluctuations is the ballistic time 

Tcl = Teol = icol/WE- 



9 



Having characterized the fluctuations, we can deter- 
mine the bandwidth = h/r^^i. A straightforward cal- 
culation leads to the result: 



L 



(26) 



where Ae — 27:%/ {mvj^) is the De-Broglie wavelength, 
and L± is the width of the ring. Using Eq.|(7|) we can 
determine the non-perturbative parametric scale: 



L 



(27) 



which up to a geometric factor equals the quantal flux 
unit. Note that in order to mix levels a relatively small 
change in the flux is needed, as implied by comparing 
Eq.® with Eq.O. 

We turn now to the analysis of the spreading in the 
time dependent case, say for linear driving. The classical 
"slowness" condition which has been mentioned in sec- 
tion 11 is simply eV <^ E where E is the kinetic energy 
of the charged particle. Upon quantization we should 
distinguish the non-perturbative regime using Eq. 
leading to 



eV > 



L 



(28) 



Disregarding the geometric prefactor, the quantity in the 
right hand side is the so called Thouless energy. We 
also should distinguish the QM adiabatic regime using 
Eq. , leading to 



(29) 



where we have assumed for simplicity a simple 2D quan- 
tum dot geometry as in Fig. 4a. 



XIV. DRIVING BY MOVING WALLS 

There is an ongoing interest [TEl IT^ in the problem 
of ID box with moving wall (also known as the infinite 
well problem with moving wall). If the wall is moving 
with constant velocity, then it is possible to transform the 
Schrodinger equation into a time-independent equation, 
and to look for the stationary states. 

We are interested in the dynamics, and therefore we 
would like to go beyond this limited scope of study. Be- 
fore we discuss the general case, it is useful to point out 
the d > 1 generalization of the above picture. We can 
easily show that for any "special deformation" which is 
executed in either constant "velocity" or "acceleration", 
we can transform the Schrodinger equation into a time- 
dependent equation. By "special deformation" we mean 
either translation or rotation or dilation, or any linear 
combination of these. The statement is manifestly trivial 



for translations and rotations (it is like going to a differ- 
ent reference frame), but it is also true for dilations. The 
ID box with moving wall is just a special case of dilation. 

It is important to realize that in case of generic defor- 
mation of chaotic box, we cannot "eliminate" the time 
dependence. Thus it is not possible to reduce the study 
of " dynamics" to a search for " stationary solutions" . 

The determination of for this system is quite obvi- 
ous but subtle lU . As one can expect naively the result 
is Ah — 27r?i/Tcoi, where TcoI is the mean time between 
collisions with the moving walls. The subtlety here is 
that we cannot interpret A^ as "bandwidth". Formally 
the correlation time of T{t) is — oo which implies infi- 
nite bandwidth. Still, some non-trivial reasoning [30ll3l| 
leads to the conclusion that the naive result (rather than 
the "formal" one) is in fact effectively correct. A straight- 
forward calculation leads to the result: 



A 



(30) 



where Ae = 2'Kfi/ {mvj^) is the Dc-Broglie wavelength, 
and Vbox is the volume of the box, and Leoi is the mean 
free path between collisions. As for the effective value 
of (Jajprt , again the details are subtle, but the naive guess 
turns out to be correct. With the proper (natural) choice 
of units for the displacement parameter x, the result is 
simply fep.t = Ae- 

The way to analyze the dynamics for box with mov- 
ing walls is outlined in 62]. The classical LRT domain 
\s V We, where = \/2E/m. Upon quantization 
we should distinguish the non-perturbative regime using 
Eq.(j21|), leading to 



V > 



(31) 



In the non-perturbative regime the dynamics has a semi- 
classical nature, and the energy spreading process has a 
resonant random- walk nature. This should be contrasted 
with the behavior in the perturbative non-adiabatic 
regime, where Fermi-golden-rule (FGR) picture applies. 

We also should distinguish the QM adiabatic regime 
using Ea. H25|) . leading to 



V < 



3/2 



(32) 



In the QM adiabatic regime the spreading is dominated 
by transitions between near- neighbor levels: This is the 
so called Landau-Zener spreading mechanism 4]. See 
also Section 20 of [l^, and the numerically related work 
in |50|. 
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XV. BROWNIAN MOTION 

Brownian motion is a paradigm for the general prob- 
lem of system that interacts with its environment. (See 
Fig. 15 and general discussion in the next section). One 
can imagine, in principle, a "zoo" of models that describe 
the interaction of a Brownian particle with its environ- 
ment. However, following Caldeira and Leggett j5^], the 
guiding philosophy is to consider "ohmic models" that 
give Brownian motion that is described by the standard 
Langevin equation in the classical limit. Four families of 
models are of particular interest: 

• Interaction with chaos. 

• Interaction with many-body bath. 

• Interaction with harmonic bath. 

• Interaction with random-matrix bath. 

Below we assume that the total Hamiltonian has the fol- 
lowing general form 

- noix,p) + n{Q,P;x) (33) 

where (x,p) are the system coordinates, and {Q,P) are 
the environmental degrees of freedom. 

Interaction with chaos provides the simplest model for 
Brownian motion 51] . See Fig. 16a for illustration of the 
model. The large Brownian particle is described by the 
canonical coordinates {x,p), while the gas particles are 
described by the canonical coordinates {Q,P). It is im- 
portant to realize that in order to have Brownian motion, 
it is not essential to consider " many particle gas" . " One 
particle gas" in enough, but the motion of the gas particle 
should be chaotic. 

The fluctuations of the environment are in fact (ac- 
cording to our definition in Section 3) the random-like 
collisions of the gas particle with the Brownian parti- 
cle. These fluctuations are like "noise". If the motion 
of the gas particle is strongly chaotic, then the power 
spectrum of these fluctuations (Fig. 5) is just like that of 
white noise. [This characterization is meaningful up to a 
cutoff frequency which is determined by the rate of the 
collisions.] 

On the other hand we have the effect of dissipation. If 
the particle is launched with a velocity i — V, then the 
rate of dissipation is proportional to as explained in 
section 9. Having dissipation implies that the Brownian 
particle experiences friction force which is proportional 
to the velocity V. This is the reason why the dissipation 
coefficient is known also as friction coefficient. 

Interaction with chaos can be regarded as the "meso- 
scopic" version of Brownian motion. Our interest in this 
set of lectures is in this type of interaction. We want 
to know whether few degree of freedoms can serve as a 
"bath". Before we further get into this discussion we 
would like to describe the "conventional" point of view 
regarding Brownian motion. The rest of this section is 
dedicated for this clarification. 



The conventional point of view regarding Brownian 
motion assumes an interaction with many body bath. 
We can consider a bath that consists of either Bosons or 
Fermions j56l l57j . The emerging models are quite com- 
plicated for analysis, and therefore, as already mentioned 
above, it is more common to adopt a phenomenological 
approach. 

Interaction with (many body) harmonic bath is not 
very natural, but yet very popular model for Brownian 
motion. In order to have "white noise" (at high temper- 
atures or in the classical limit) we should make a special 
assumption regarding the frequency distribution of the 
bath oscillators. This is known in the literature as the 
"ohmic choice". [The characterization of the noise as 
" white" is valid up to some cutoff frequency. The latter 
is determined by the specific choice of the frequency dis- 
tribution.] Also here, as in the case of interaction with 
chaos, we have fluctuation-dissipation theorem that im- 
plies "ohmic" dissipation rate (proportional to V^). 

There is still some freedom left in the definition of the 
interaction with the harmonic bath. This leads to the 
introduction of the Diffusion-Localization-Dissipation 
(DLD) model [Ulillll]. This model gives in the clas- 
sical limit Brownian motion which is described by the 
standard Langevin equation (white noise -t- ohmic dis- 
sipation}. The familiar Zwanzig-Caldeira-Leggett (ZCL) 
model |j53 can be regarded as a special limit of the DLD 
model. The physics of the ZCL and of the DLD model 
is illustrated in Fig. 16b and Fig. 16c respectively, and the 
model Hamiltonians can be visualized by the drawings of 
Fig. 17. The ZCL model describes a motion under the in- 
fluence of a fluctuating homogeneous field of force which 
is induced by the environmental degrees of freedom. In 
case of the DLD model the induced fiuctuating field is 
further characterized by a finite correlation distance. 

For completeness we note that random-matrix model- 
ing of the environment, in the regime where it has been 
solved (53, leads to the same results as those obtained 
for the DLD model. 



XVI. SYSTEM INTERACTING WITH 
ENVIRONMENT 

The general problem of system that interacts with its 
environment is of great importance in many fields of 
physics. The basic ingredients of this interaction are il- 
lustrated in Fig. 15. On the one hand we have the effect 
of dissipation, meaning that energy is lost by the "sys- 
tem" (Brownian particle) and is absorbed by the "en- 
vironment" (gas particles). On the other hand the en- 
vironment induces fiuctuations that acts like "noise" on 
the system. The "noise" has two significant effects: One 
is to pump "thermal" energy into the system, and the 
other is to spoil quantum coherence. The latter effect is 
called decoherence. 

In case of bounded system, in the absence of external 
time dependent fields, the interplay between "noise" and 
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"dissipation" leads eventually to "thermalization" . One 
may say that in the thermal state the effect of dissipa- 
tion is balanced by the energy which is pumped by the 
noise. Thus, both classically and quantum mechanically 
we have to distinguish between a " damping" scenario and 
an "equilibrium" situation. The thermalization process 
is traditionally described as " irreversible" . On the other 
hand we have the issue of " recurrences" . We discuss the 
latter issue in Section 20. 

A systematic approach for the study of the dynamics 
of a " system" , taking into account the influence of its en- 
vironment, has been formulated by Feynman and Vernon 
^9J. The state of the system is represented by the prob- 
ability matrix p{x',x"). It is assumed that initially the 
" system" is prepared is some arbitrary state. Its state at 
a later time is obtained by a propagator K{x' , x"\x'q, Xq) 
which acts on the initial preparation. The calculation 
of this propagator involves a double path integral over 
all the possible trajectories XA{t) and xsit) that connect 
{x'q,Xq) with {x',x"). This double path integral incor- 
porates the effect of the environment via an "influence 
functional" which is defined as follows: 

F[xa,xb] - {U[xb]'^o\U[xa]-9o) (34) 

Here 5*0 is the initial state of the environment. If the 
environment is in "mixed" state, typically a canonical 
(thermal) state, then the influence functional should be 
averaged accordingly. 

The absolute value of the influence functional can 
be re-interpreted as arising from the interaction with 
a c- number noisy field (with no back reaction). The 
"phase" of the influence functional can be regarded 
as representing the effect of "friction" (back reaction). 
Thus there is one to one correspondence between the 
Feynman- Vernon formalism, and the corresponding clas- 
sical Langevin approach. Note however that the distinc- 
tion between " noise" and " friction" is a matter of " taste" . 
Some people regard this distinction meaningless. 

It should be realized that the calculation of the influ- 
ence functional for a given environment takes us back to 
the more restricted problem of considering a " driven sys- 
tem". The influence functional F[xa,xb] is nothing but 
the survival amplitude for a driving scheme that involves 
"time reversal" fEa. lfT^ l. 

XVII. ENTANGLEMENT, DECOHERENCE 
AND IRREVERSIBILITY 

The definition of decoherence is not a trivial matter 
conceptually. There are several equivalent ways to think 
about decoherence. The most "robust" approach is to de- 
fine decoherence as the irreversible entanglement of the 
system with the environment: Let us describe the the 
state of the system using the probability matrix p{x' , x"). 
If the system is prepared in pure state then trace(/9^) = 1. 
Due to the interaction with the environment the sys- 
tem gets entangled with the environment, and as a re- 



sult we will have trace(p^) < 1. If the "environment" 
consists of only "one spin", then we expect trace(p^) to 
have "ups" and "downs", and from time to time to have 
trace(p^) ~ 1. In such case we cannot say that the entan- 
glement process is " irreversible" . But if the environment 
consists of many degrees of freedom, as in the case of in- 
teraction with "bath", then the loss of "purity" becomes 
irreversible, and we regard it as a "decoherence process". 

To be more specific let us consider the prototype ex- 
ample of interference in Aharonov-Bohm ring geometry. 
The particle can go from the input lead to the output 
lead by traveling via either arms of the ring. This leads 
to interference, which can be tested by measuring the 
dependence of the transmission on the magnetic flux via 
the ring. Consider now the situation where there is a 
spin degree of freedom in one arm [23| . The particle can 
cause a spin flip if it travels via this arm. In such case 
interference is lost completely. However, this entangle- 
ment process is completely reversible. We can undo the 
entanglement simply by letting the particle interact with 
the spin twice the time. Therefore, according to our re- 
strictive definition, this is not a real decoherence process. 

Consider now the situation where a particle gets en- 
tangled with bath degrees of freedom. If the bath is in- 
finite, then the entanglement process is irreversible, and 
therefore it constitutes, according to our definition, a de- 
coherence process. 

At first sight it seems that for having irreversibility 
one needs " infinity" . This point of view is emphasized 
in Ref. [6Clt | : Irreversibility can be achieved by having the 
infinity of the bath (infinitely many oscillators), or of 
space (a lead that extends up to infinity). In this set of 
lectures we emphasize a third possibility: Having irre- 
versibility due to the interaction with chaos. Thus we do 
not need " infinity" in order to have " irreversibility" . 



XVIII. INTERPRETATION OF 
DECOHERENCE AS A DEPHASING PROCESS 

"Dephasing" is used as a synonym for "decoherence" 
whenever a semiclassical point of view is adopted. De- 
termining the dephasing (decoherence) time for a par- 
ticle {x,p) that interacts with an environmental degrees 
of freedom (Q, P) is a central theme in quantum physics. 
In the absence of such interaction the x motion is coher- 
ent, and interference should be taken into account. This 
means, from semiclassical point of view, that at least two 
trajectories x{t) = x^^t) and x{t) = Xb{t) have a lead- 
ing contribution to the probability to travel, say, from 
x{0) to x(t), as in the prototype example of the two slit 
experiment. 

In the semiclassical framework the probability to travel 
from one point to some other point is given by an expres- 
sion that has the schematic form 

1 f .S[xa]- S[xb]\ 
2_^F[xA,XB]eTip It I (35) 

A,B \ n- J 



12 



where S[x] is the classical action, and F[xa,xb] is the 
influence functional. Each pair of trajectories is a "sta- 
tionary point" of the Feynman- Vernon double path in- 
tegral. The "diagonal terms" are the so-called classical 
contribution, while the "off-diagonal terms" are the in- 
terference contribution. It should be kept is mind that 
the validity of the semiclassical approach is a subtle issue 

ill- 

The off-diagonal interference contribution is sup- 
pressed due to the interaction with the environment if 
|F[a:A, ics]! 1. Therefore |F[a;A, a;_B]| is called the 
" dephasing factor" . From the definition of the influ- 
ence functional it is clear that it reflects the proba- 
bility to "leave a trace" in the environment. Having 
|F[a;A,a;B]| = means that a different "trace" is left 
in the environment, depending on whether the particle 
goes via the trajectory XAit) or via the trajectory XB{t)- 
In such case one can regard the interaction with the en- 
vironment as a "measurement" process. In case of the 
DLD model (see Section 15) this "trace" can be further 
interpreted as leaving an excitation along the way. For 
critical discussion of this point see Appendix C of [H3 |. 
In the more general case the notion of " leaving a trace" 
does not have a simple meaning. All we can say is that 
decoherence means that the environment is left in differ- 
ent (orthogonal) states depending on the trajectory that 
is taken by the particle. 

The law of "action and reaction" holds also in the 
world of decoherence studies. Feynman and Vernon have 
realized that the dephasing factor can be re-interpreted 
as representing the effect of a c-number noise source (see 
section 16). From this point of view the decoherence 
is due to the "scrambling" of the relative phase by this 
noise. Hence the reason for using the term "dephasing" 
as a synonym for " decoherence" . The analysis of dephas- 
ing using this latter point of view can be found in |54j . 
See also At high temperatures it is possible to use 
a Markovian master equation approach (dynamical semi- 
groups) in order to obtain the (reduced) evolution of the 
Brownian particle. The Markovian master equation ap- 
proach is described in other lectures of this school. The 
master equation in case of the DLD model can be found 
in Section 3 of |54|. Similar, but not identical master 
equations are obtained in case of interaction with many 
body bath 



XIX. DETERMINATION OF THE DEPHASING 
TIME 



In the above described semiclassical framework, the 
problem of dephasing reduces to the more restricted 
problem of studying the dynamics of a time dependent 
Hamiltonian H{Q,P;x{t)). Moreover, we see that the 
Feynman- Vernon dephasing factor is just the absolute 
value of the fidelity amplitude F{t) that corresponds to 
Ea. (|15|l . Note however that here we use a more gen- 
eral notion of fidelity: The restricted definition of fidelity 



fEa. (|17|l ) is formally obtained if Xa{t) and Xb(t) are 
" rectangular pulses" . 

The dephasing time is defined as the time that it 
takes for \F{t)\ to drop significantly from \F{t)\ ^ 1 to 
some very small value \F{t)\ ^ 1. Let us concentrate on 
the Brownian motion model of Fig. 16a. If the motion of 
the Brownian particle is characterized by a velocity V, 
then we have to distinguish between the following possi- 
bilities: Having very small "adiabatic" velocities; Having 
intermediate velocities that allow LRT treatment; And 
having non-perturbative velocities. In the latter case a 
semiclassical picture can be justified. 

The detailed analysis of the problem can be found in 
[6^ . Here we just quote the final results. In the semi- 
classical regime 



ice 

Ve 



(36) 



where L^^i is the mean free path between collisions with 
the Brownian particle. This is the expected naive result. 
It means that one collision with the Brownian particle is 
enough in order to "measure" its trajectory. The other 
extreme case in having extremely small adiabatic veloci- 
ties. To the extend that we can trust adiabaticity there 
is no dephasing at all: The gas particle simply "renor- 
malize" the bare potential, which is in fact the Born- 
Oppenheimer picture. Of course, if we take into account 
corrections to the adiabatic picture, then we get a finite 
dephasing time. In the LRT regime of velocities we can 
estimate the dephasing time as 



(37) 



Both results have re-interpretation within the framework 
of an effective DLD/ZCL model. See |63 for details. 



XX. RECURRENCES 

Consider ice-cube inside a cup of hot water. After some 
time it melts and disappears. But if we wait long enough 
(without time reversal) we have some probability to see 
the ice-cube re-emerging due to recurrences. The issue 
of recurrences becomes relevant whenever we consider a 
closed (un-driven) system. In other words, whenever we 
do not try to control its evolution from the outside. 

There are recurrences both in classical and quantal 
physics. In the latter case the tendency for recurrences 
is stronger due to the quasi-periodic nature of the dy- 
namics. However, if the time scale for recurrences is long 
enough with respect to other relevant time scales, then 
we can practically ignore these recurrences. Actually it is 
useful to regard these recurrences as "fluctuations", and 
to take the standpoint that our interest is only in some 
"average" or "likely" scenario. 

The thermalization process of the particle-environment 
system is traditionally described as " irreversible" . In- 
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deed, if the bath is infinite, then also the time for recur- 
rences of the particle-bath system becomes infinite. On 
the other hand, if the bath is finite, then we have to con- 
sider the recurrences of the particle-bath system. These 
recurrences can lead back to the initial un-entangled 
state. 

In practice "recurrences" do not constitute a danger 
for " irreversibility" . The time to get un-entangled by 
recurrences is extremely large (typically larger than the 
age of the universe). Assuming a chaotic environment, 
and ignoring issues of level statistics, the time scale for 
recurrences is at least the Heisenberg time (inverse of the 
mean level spacing) of the combined particle-environment 
system. It scales like where do and d are the 

number of degrees of freedom of the particle and the en- 
vironment respectively. 

It goes without saying that the above issue of recur- 
rences becomes irrelevant if the x motion is treated clas- 
sically. There is however a twist to this latter statement 
in the case where the time variation of x is strictly pe- 



riodic. This is due to dynamical localization effect 
Note however that dynamical localization is a very frag- 
ile effect: Even in case that it is found, it turns out that it 
manifests itself only after a time that scales like fi~^'^^'^'^\ 
which is much larger than the Heisenberg time of the en- 
vironment pTl |. 
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x(t) 

FIG. 1: The prototype piston model. A gas particle is mov- 
ing chaotically inside a cylinder. The driving is achieved by 
moving a wall element ("piston"). 



ID Box 



2D Box 




FIG. 2: The ID version of the piston model (upper panel). 
The gas particle is moving inside an "infinite well" . Its motion 
is not chaotic. In order to have chaotic motion we should 
consider at least a 2D box, for example a stadium shaped 
billiard system (lower panel). 




FIG. 3: Other versions of the "piston model". Two examples 
of "Sinai billiards" arc illustrated. Note that in case of the 
right panel the displacement of the wall element manifestly 
does not involve a change of volume. This feature simplifies 
the analysis. 




FIG. 4: The same "Sinai billiards" as in the previous figure. 
Here the driving is achieved by changing the perpendicular 
magnetic field. In case of the left panel, where the box has 
a simple "quantum dot" geometry, the magnetic field is as- 
sumed to be homogeneous. In case of the right panel, where 
the box has a aharonov-Bohm ring topology, the magnetic 
flux is assumed to be concentrated in the hole. 




FIG. 5: Chaotic trajectories inside a Sinai billiard. In both 
examples the motion is completely chaotic. This mean ex- 
ponential sensitivity to any small change in the initial condi- 
tions. This sensitivity can be characterized by the "Lyapunov 
exponent" . In the right illustration the motion is chaotic, but 
the chaos is weaker, which means smaller Lyapunov exponent. 
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FIG. 6: The power spectrum of the motion for the two ex- 
amples of the previous figure. It is the power spectrum C{u}) 
of the fluctuating quantity J-{t) = —dTL/dx. The latter can 
be described as train of impulses (spikes) due to collisions 
with the walls. For strongly chaotic motion (left panel) the 
power spectrum of J-{t) is like that of white noise. In the 
right panel the bouncing frequency is quite pronounced, and 
there is also a "diabatic" peak around = 0. In both cases, 
the motion is characterized by a continuous power spectrum, 
which constitutes an indication for the chaotic nature of the 
motion. 
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classical 
h=0.030 
h=0.015 




FIG. 7: The band profile (27r;i/A) • |B„™p versus uj = 
{En—Em)/fi- is compared with C{ijj)- See text for further 
explanations. The calculation is done for the 2DW model of 
Eq.ll^J. The inset is an image of a piece of the B matrix. 
Taken from Ref.H. 
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FIG. 8: Upper panel: Schematic illustrations of the set of 
energies En{x) which are obtained via diagonalization of a 
parameter dependent Hamiltonian. The thick solid line indi- 
cates the n-range where 50% of the P{n\m) probability is con- 
centrated (m is fixed). The representative values xi, X2 and 
xq correspond to the standard perturbative regime, the core- 
tail (extended perturbative) regime, and the non-perturbative 
regime respectively. The corresponding LDOS structures are 
illustrated (grey shading) in the three plots of the lower panel. 
The semiclassical approximation (lines) is presented for sake 
of comparison. 
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FIG. 9: Various driving schemes: Rectangular pulse of du- 
ration T and amplitude A; Triangular pulse which is further 
characterized by finite driving "velocity" V — \x\ = 2A/T; 
Periodic driving with frequency Q = 2-k/T. 
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FIG. 11: The dependence of the friction coefficient on the 
driving frequency, using LRT. " GN" stands for generic defor- 
mation of a stadium shaped billiard, while "DI" stands for 
special deformation (dilation). [In the latter case the fric- 
tion coefficient vanishes in the low frequency limit.] In both 
cases the agreement between the classical (solid line) and the 
quantum-mechanical (dashed line) calculation is remarkable. 
Taken from Ref.^. 




FIG. 10: Energy spreading as a function of time for the 2DW 
model: (a) classical; (b) quantum mechanical; (c) An effective 
Wigner model. The energy in these simulations is -E ~ 3, and 
Sx = 0.2123. The light dashed lines in (a) and (c), that have 
slopes 1 and 1/2 respectively, are drawn to guide the eye. In 
(c) different lines correspond to different values of h as in (b), 
and additional curves {h = 0.009, 0.005) have been added. 
Taken from Ref.ilJl. 
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FIG. 12: The dependence of the mesoscopic conductance on 
the driving frequency. The calculation is done for a Sinai 
billiard shaped quantum dot, using LRT. The result can be 
regarded as a mesoscopic version of Drude formula. The spe- 
cific geometry of the system is reflected in the structure of the 
response curve. The inset is log-log plot. Taken from Ref . . 
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FIG. 13: The response of a quantum mechanical system is 
displayed as a function of A and Q.. The simulations are done 
for Wigner model with periodic driving. The units are chosen 
such that A = 0.5 and fi = 1 and ct = 1. Upper panel: Plots 
of D-e/A^ versus O/oioi for few values of A. For small lo the 
plots coincide as expected from LRT. As A becomes larger 
the deviations from LRT scaling become more pronounced, 
and we get response also for Q. > Wci. Lower panel: Plots 
of De/Dq versus Ajyfh for few values of O/i^ci. LRT implies 
-De /Do = 1 for ^ju^i < 1 and De/Do = for Q/lOoi > 1. The 
purpose of the horizontal scaling is to demonstrate that A-p^t 
rather than Ac is responsible for the deviation from this LRT 
expectation. Taken from R.ef . [TH . 
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FIG. 14: Upper diagram: The various V regimes in the the- 
ory of driven systems for linear driving x{t) = Vt. Lower 
diagram: The various (fi, A) regimes for periodic driving 
x{t) — Asm{flt). We use the notations cJci = 27r /tci and 
Ac — Sxc and Aprt = Sxp^f 
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FIG. 15: Block diagram illustrating the interaction between 
system (x) and environmental (Q) degrees of freedom (DoF). 
See discussion in the text. 
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FIG. 16: (a) The Brownian motion is induced due to the 
interaction with chaotic degrees of freedom, (b) The Brown- 
ian particle in the ZCL model experiences a fluctuating ho- 
mogeneous field of force, (c) In case of the DLD model the 
fluctuating field is farther characterized by a finite correlation 
distance. In (b) and (c) the background image is a "snapshot" 
of the fluctuating environment. Namely, the gray levels cor- 
respond to the "height" of an instantaneous potential which 
is experienced by the Brownian particle. 



FIG. 17: Schematic illustration of the ZCL model (upper 
panel) and the DLD model (lower panel). The Hamiltonian 
of these "spring systems" is literally the ZCL model and the 
DLD model respectively. In the latter case the height of the 
masses should be interpreted as the "height" of the potential 
which is experienced by the particle. 




FIG. 18: Schematic illustration of an interference experiment 
using a semiclassical point of view. The Brownian particle 
can take either the x = XA{t) trajectory, or the x = XB(t) 
trajectory as in a two slit experiment. 



